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abstract 

We investigate a symmetry broken dimer constrained to move in a particular direction when in 
contact with a heat-bath at a constant temperature. The dimer is not driven by any external force. 
The system gains kinetic energy from the heat-bath. The symmetry broken system can use this 
energy in directed transport. At the hard core collision limit between the particles of the dimer, we 
show by exact analytic calculations and complementary numerical results that the dimer undergoes 
steady directed transport. Our observation, being consistent with the second law of thermodynamics, 
detailed balance etc leads to new physical understanding to which much attention has not been paid. 
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Introduction 

The Brownian Ratchet (BR) paradigm provides founda- 
tion to the understanding of symmetry-broken driven di- 
rected transport in contact with a heat-bath 0-11]. BR 
paradigm is based on the concepts of Smoluchowski- 
Feynman ratchet that illustrates a typical case where ex- 
tracting energy from the random motions of a heat-bath 
alone is impossible due to violation of the second law 
of thermodynamics Jy\. In general, BR paradigm claims 
that thermal noise assists driven transport in the pres- 
ence of a repeated spatial asymmetry. Attempts have 
been made in understanding biological transports in the 
light of BR paradigm by invoking numerous variants of 
the basic BR model in order to match the efficiency of the 
actual biological systems [6-14]. However, the efficiency 
matching is quite poor till date. There are many other 
areas, like flow through micro-channels, nanopores, opti- 
cal lattices etc., where the BR paradigm finds its appli- 
cation to the noise assisted driven transport phenomena 
0, Il5l - [l7j ]. There also exists a variant of conventional 
BR, called an Intrinsic Ratchet (IR) model, where the 
internal degree of freedom is not coupled to an exter- 
nal global symmetry breaking field. Interestingly, such 
IR systems with microscopic broken symmetry are ca- 
pable of showing macroscopic directed transport in the 
presence of a drive that disturbs the equilibrium of the 
system intermittently lil 2pj] . 

It is generally believed that, the Smoluchowski- 
Feynman ratchet model denies any filtering of random 
motions of a heat bath under equilibrium conditions. 
Moreover, when the dynamics of an interacting system 
of particles is cast in the (over-damped) Langevin equa- 
tion (LE) form 2ll~l23j. the position coordinate becomes 
the relevant variable. The demand of a stationary proba- 
bility distribution of the variable of this general LE, as a 
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condition for the equilibrium of the system with the heat- 
bath, would (arguably) preclude any directed transport 
of the system. Because, any average directional drift of 
the system might mean a drift in the stationary probabil- 
ity distribution of the position coordinates which will be 
in conflict with the demand of stationary distribution of 
position variables. But, as one can reasonably guess, (i) 
the second law of thermodynamics would not forbid such 
a constant drift so long as no energy is being extracted 
from the system, (ii) the detailed balance on a uniformly 
moving frame is maintained due to homogeneity of space 
which makes interactions between particles functions of 
their distance of separation and (iii) nonexistence of time 
reversal invariance for a moving system is not inconsis- 
tent with second law of thermodynamics (Loschmidt's 
paradox). To our knowledge, possibility of such directed 
transport of a mechanical system in equilibrium with a 
heat-bath (needless to say in the absence of any external 
driving) has not been systematically explored so far. In 
the present paper, we investigate one generic model for 
directed transport in equilibrium with a heat bath. 

We consider a microscopic system with internally bro- 
ken symmetry like the IR, but, unlike IR its not driven for 
directed transport. Rather, we investigate a set of com- 
pletely different issues to understand what happens when 
such a symmetry-broken microscopic object is brought in 
contact with a heat-bath and the locally broken symme- 
try is not allowed to relax. If the origin of the broken 
internal symmetry of the system is consistent with condi- 
tions of equilibrium with the heat-bath, the system would 
definitely have a preferred direction of motion and must 
also equilibrate at large times. If it does not equilibrate, 
that is a violation of second law of thermodynamics in it- 
self, because, under such ever-persisting non-equilibrium 
situation one may extract energy from a heat-bath for 
ever without external driving. In the following, we first 
identify a general situation in which such a system would 
show directed transport. We show exact analytic results 
for elastic hardcore collisions and provide numerical sup- 
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port to it. We identify the damping or dissipation as the 
key physical ingredient required for symmetry breaking 
directed motion in equilibrium. We conclude with a dis- 
cussion of our observations in relation with second law, 
detailed balance, entropy production etc. 
Model 

The dynamics for our model system is governed by the 
set of equations 
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The mass of the particles at x\ and x 2 [x\ > x 2 ) are 
unity and damping constants are (1 — ft) and 1 respec- 
tively. The damping can be different on different par- 
ticles depending upon size, shape etc. of the particles, 
where by particle we do not mean dimensionless point 
objects. The random forces on these particles rji and r) 2 
are Gaussian white noise with the moments (T]i(t)) = 
and (r]i(ti)rij(t 2 )) = Sijd(ti—t 2 ). The particles are bound 
by a spring force of force constant a. In the last term 
$(xi —x 2 ) is the collision potential taking care of the ex- 
cluded volume interaction keeping x\ > x 2 for ever. The 
strength of the thermal noises are explicitly mentioned 
in accordance with the fluctuation-dissipation relation so 
that the particles see the same temperature T of the heat- 
bath where ks is the Boltzmann constant. Similar mod- 
els with broken symmetry due to different damping have 
been considered in ref.[2J, |25|, but, systems are driven 
there. 

Let us think of a situation where the collision potential 
$>(x% — x 2 ) is so steep that the actual collision time (the 
time in which particles are in contact with each other) is 
very very small compared to the slow time scale over 
which the over-damped dynamics of the system takes 
place. This situation allows us to separate the dynamics 
into two parts. The inertial part only takes care of the 
collision forces and the over-damped part takes care of 
the rest. The inertial part on the center of mass (CM) 
X = (xi + x 2 )/2 and internal coordinate Z = x\ — x 2 
would look like 
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The inertial part does not contain other forces because 
of the following reasons. The damping is particularly 
defined at a larger time scale and there is every reason 
to take the stochastic force averaging out over several 
such intermittent collisions because they are not corre- 
lated. In the following where we would be introducing 



the corrections occurring from the collisions to the av- 
erage velocity of the CM and Z, we would actually be 
introducing the average contributions over several such 
collisions. Moreover, this is a plausible hypothesis be- 
cause we can actually keep a, in the present model, ar- 
bitrarily small which would in turn make average time 
gap between collisions large and that can actually make 
the stochastic contribution at collisions uncorrelated. So, 
based on the above logic of considering time scale corre- 
sponding to the harmonic interaction large Eq.4 clearly 
indicates that there will be no contribution to the CM 
velocity from the collision whereas Eq.3 indicates a con- 
tribution to the relative velocity of the particles coming 
from the collision which would be just in the opposite di- 
rection to that coming from the harmonic part. For the 
system to remain in equilibrium the < Z > must vanish 
and that would happen due to cancellation of contribu- 
tions coming to it from the harmonic interaction and this 
collision part. In what follows we would take care of this 
short time scale dynamics through the boundary condi- 
tions imposed to the large time scale dynamics. Its im- 
portant to note that Eq.l and 2 with all the forces present 
simultaneously does not represent our model. Our model 
particularly relies on the non-simultaneous acting of the 
collision and damping like slow forces. This separation of 
dynamics is crucial to get the expected results because it 
will be shown in the following that the effective symmetry 
breaking of the slow harmonic interaction by the damp- 
ing keeping the collision symmetry intact is the basis of 
getting the directed motion of the model. It is also impor- 
tant to note, at this stage, that this somewhat idealized 
model we are going to analyze is entirely because of the 
simplicity. Based on this model, we would understand a 
more realistic situation taking into account inelastic col- 
lisions of the particles where we can actually have such 
conditions of idealized hardcore collisions relaxed. 

The equilibrium of different particles (in shape/size 
etc.) with the same heat-bath ensures broken symme- 
try and this is a key aspect of the construction of the 
present model. Considering the over-damped (slow) part 
of the dynamics to cast the model into a standard cou- 
pled Langevin equation 2l|, |22j and moving on to the 



internal coordinate and CM, we rewrite the above set of 
equations in a more convenient form as 
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where U(Z) = %Z 2 . The modified random forces £s 
are functions of 77s and are of zero average. Explic- 
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itly, expressions of the noises are £z = 
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The second moments of the random forces are given as 
(£z(*i)fc(*a)) = 2Tk B fE$S(h - t 2 ), (ex(ti)6c(* a )) = 
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£2)- The steady state probability distribution for the 
Z, when f7(i?) is well behaved, can be readily given as 
P{Z) = N cxp(-U(Z)/k B T) (N is normalization con- 
stant) with a sharp fall to zero value at Z = Zi which 
is the boundary condition imposed in view of the steep 
excluded volume interaction (fast collision process) . This 
probability distribution can be used in general to write 
down an expression for the average velocity of the CM as 
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Important to note that, in the above integral we would 
not take into account the delta function contribution that 
comes near Z = Zi. This delta function contribution rep- 
resents the contribution of the collision potential to the 
velocity of CM which we know from Eq.4 to be identically 
equal to zero. The only contribution that would come to 
the velocity of the CM is from the harmonic part of the 
interaction. Remember that, its the interactions over the 
span of the harmonic potential only which contributes 
to the development of the velocity of the CM and dur- 
ing the collisions between the particles constituting the 
dimer that velocity remains conserved. But, for < Z > 
we must take into account the delta function contribution 
to the integral at Z = Zi which would make < Z >= al- 
though < X >t^ 0. We would take into account the delta 
function contribution for the average Z because this part 
representing the collision indeed contributes to the rel- 
ative velocity of the particles of the dimer with respect 
to each other. Its interesting to note that, consideration 
of a very steep collision potential amounts to abruptly 
bringing the P(Z) from its maximum value to zero over 
a very short range of Z near Z = Z\ which is at the ori- 
gin of consideration of the delta function contribution. 
Also note that, since only the harmonic interaction is 
contributing to the velocity of the CM, < X > is propor- 
tional to < Z >. The reason the harmonic interaction 
is contributing to the average CM velocity is that its ef- 
fectively asymmetric and the symmetries of the collision 
interaction prevents it from having any contribution to 
the CM velocity. 
Results 

From the over-damped dynamics one obtains an equi- 

-c,Z 2 

librium distribution of Z as P{Z) = Ne 2k B T within 
the range < Z < oo (here Zi = and at this point 
P(Z) = is the boundary condition). If the system in- 
deed attains this steady state distribution at large times 
from any initial conditions, the steady average velocity 
of the CM will be given by 
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FIG. 1. Average velocity V of CM plotted against 
force constant a and damping constant (3: (a) Average 
velocity vs a at ft — 0.5 and (b) Average velocity vs /3 at 
a = 1.0. 



which is an exact result. 

A direct numerical simulation of Eq.5 and 6 with this 
elastic hardcore collision is actually possible and reveals 
numbers in excellent agreement with the exact analytic 
results. In Fig. 1(a), we have plotted the average veloc- 
ity of CM against the parameter a and in Fig. 1(b) we 
have plotted the same against /3. While doing the nu- 
merical simulation, we have implemented the repulsion 
between particles with the following rules. Since, for par- 
ticles of same mass (as has been considered here), there 
is an exchange of velocity at elastic hard-core collisions, 
we simply exchange the position of the particles follow- 
ing any overshooting across each other during the sim- 
ulation. This simple implementation of boundary con- 
ditions, while sequentially updating particle positions, 
takes care of the perfect hard-core collision preserving 
the momentum at collisions. Similar agreements have 
also been seen at the perfect inelastic collision limit where 
we put both the particles at the CM following any over- 
shooting and allow them to move in their respective sides 
only in the next time step if no overshooting is happening 
again. The perfect inelastic collision also showing simi- 
lar results due to the fact that it also does not mess up 
with the CM momentum conservation at collisions. The 
case of this inelastic collision is the key to understand 
the physically realizable situation for such motions. We 
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FIG. 2. Comparison of numerical and analytical probability 
distribution of Z: Distribution of Z for a = 1.0 at = 0.5. 

will discuss about it at length in the following, but, here 
we would like to mention that under the scope of our 
present model, where the temperature of the bath is kept 
constant despite implementing inelastic collision between 
particles in the simulation, we are considering that the 
particles are relaxing back quickly by giving back the en- 
ergy to bath following an inelastic collision. The excellent 
agreement between the theory and the simulation indi- 
cates that the system has indeed attained its equilibrium 
distribution in Z which clearly comes out to be the case in 
Fig. 2 where numerical and analytical distributions of Z 
are compared. We have checked this comparison of the- 
ory and simulation for a range of a and (3 and those are 
quite consistent. The excellent matching of probability 
distributions indicates that corresponding Fokkcr-Planck 
equation (FPE) works fine for a U(Z) having discontinu- 
ity at the boundary. In Fig. 3 we have plotted few trajec- 
tories of the moving dimer for the illustration of the av- 
erage directed motion of the system at various parameter 
values. In Fig. 5 we have shown a comparison of < Z > 
against a as obtained from theory and numerics. In all 
these simulations the ksT = 0.00005, the time step is 
At = 0.0001, a single run is over 10 7 such time steps and 
the averaging is done over 10 such runs starting from the 
same initial conditions for each set of parameters. Note 
that, the characteristic time scale for this system being 
ft/a is of the order of 0.01. 

Consider the fact that, the dimer is maintaining a uni- 
form velocity in the presence of damping. So, what bal- 
ances this damping to help the system maintain its veloc- 
ity? Eq.6 depicts the over-damped dynamics of the CM 
in the presence of two forcing terms - the one originating 
from the harmonic interaction between the particles and 
the other one being the random force due to the bath de- 
grees of freedom and that is of zero average. Obviously, 
there is no contribution to the average velocity of the CM 
from the random force part because it has a zero average. 
Also remember that, in the case of hardcore collisions the 
contribution from the collision forces to the average CM 
velocity is zero (as we have already shown). To undcr- 




400 

t 



FIG. 3. Trajectories of the dimer for hard-core collisions: 
/3 = 0.5 and a — 0.5, 1.0 and 4.0 from top to bottom. 

stand the effective force working against the damping 
refer to Fig. 4a where we have plotted the effective spring 
potentials as seen by the particles at x\ (right) and xi 
(left). To make the picture more revealing we have plot- 
ted the potentials in such a way that the inward direction 
of the resulting forces are evident. The effective spring 
constant for the particle at x\ being (1 — j3)~ l times big- 
ger than that of the particle at x 2 , there is an effective 
force on the composite system from x\ towards x^. This 
effective force would on average continuously balance the 
damping of the CM velocity. Interesting to note that, 
the velocity picked up by the system is in the direction 
of this effective force which must have happened at the 
initial transients when the average uniform velocity was 
being reached for which the damping (which is propor- 
tional to velocity) on average is the same as the effective 
force. 

Refer to the schematic diagram of Fig. 4b to understand 
why there is no velocity of the CM when there is a soft- 
core conservative collision (collision time is of the order of 
the characteristic time scale of over-damped dynamics) 
between the particles. Here, we have considered a col- 
lision potential of the simple form $(Z) = (a — bZ) for 
< Z < e where e is the width of the collision potential. 
a and b are suitably large numbers (depending on fc^T) 
to make sure there is no passage of particles through each 
other. Moreover, we consider that the damping due to 
the bath works equally when particles actually collide 
with each other as when they are on flight. In such a 
situation, since, the attractive (harmonic) and repulsive 
(collision) effective forces between particles are similarly 
renormalizcd by their respective damping constants, in- 
dividually on average there will be a force balance for 
each of the particles. Its not difficult to understand how 
the gain of the CM velocity in one direction is counterbal- 
anced in such a situation when the collision potential is 
conservative. The velocity of approach and the velocity 
of the separation of the particles just at the point of col- 
liding would be the same and the velocity of separation 
would now counterbalance the gain in CM velocity that 
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FIG. 4. Schematic diagram of effective potentials. 



P = 0.5 




FIG. 5. Plot of < Z > against a. 

happened during the approach till the next approach is 
initiated as the particles are gone far apart. Interesting 
to note that, in this case of soft conservative collisions 
in the presence of the symmetry breaking damping, the 
velocity of separation is achieved by the velocity reversal 
and not by the velocity exchange. As a result the velocity 
of separation would contribute as much to the velocity of 
CM as the velocity of approach does, but, it would do 
that in just the opposite direction. At this point, con- 
sider what happens if the collision is perfectly inelastic. 
The velocity of separation being zero there is nothing to 
counterbalance the CM velocity gained by the help of 
velocity of approach. So, under the conditions that the 
energy absorbed by the particles at inelastic collision are 
quickly enough given back to the bath and no energy is 
lost from the composite system of the bath and the dimcr 
the system would show directed transport even for a soft- 
collision between the particles. Indeed, we have seen the 
hard-core perfectly inelastic collision to produce exactly 
the same results as the hardcore elastic collision in our 
numerical simulation. Let us discuss some subtleties of 
this inelastic collision case in the following paragraph a 
little more elaborately. 

Note that, for soft collision potentials there may come 
an explicit contribution of the collision potential in Eq.5 
and 6 because now the resulting collision forces can be a 



part of the over-damped dynamics. If we think in terms 
of the delta function contribution to the integral in Eq.7, 
then, those are now there for the velocity of the CM and 
that of the Z making both vanish. At the origin of the 
collision contribution to the average velocity of the CM 
is the effective symmetry breaking of the collision force 
by the unequal damping exactly the same way as that 
happened for the harmonic interaction. There is an way 
out though, and that involves the possibility of having 
a different damping during collisions (when particles are 
in touch with each other). Since the environment when 
the particles are at collision is different from that when 
they are in flight the damping can actually be different 
at these two times. Furthermore, there can be inelastic 
collisions which is absolutely of a different origin than the 
fluid friction and which can make the overall damping the 
system of particles experience during collisions become 
different from that when they are in flight. One has to 
keep in mind that if there is inelastic collision the energy 
ultimately has to be given up by the system to the bath 
and the system cannot behave as an infinite heat sink. 
Thus, the system has to equilibrate with the bath at large 
times but this equilibration time should now be different 
depending upon the relaxation time of the internal vibra- 
tion modes of the particles etc. Most important to note 
that, under such conditions, the fluctuation dissipation 
must hold good when the system is in equilibrium. But, 
more plausibly, it would hold good at a larger time scale 
with an average damping which is an average over the 
different damping the particles are actually undergoing 
while on flight and at collision. This is an equilibrium 
situation where there are more long lived fluctuations in- 
volved then usual due to the role played by other internal 
degrees of freedom of the system equilibrating with the 
bath. This kind of a situation, to my knowledge, is not 
quite well studied so far. This is a situation which again 
comprises of two time scales within the average time scale 
over which the FDT is defined and one can see a varia- 
tion of the damping at these two time scales. Having two 
damping, it opens up the possibility, in a similar way as 
discussed above, of having persistent broken symmetry 
in the system in equilibrium and, hence, directed motion 
for soft collisions as well. Probably this is the scenario 
more easier to probe experimentally if one can ensure the 
isolation of the system and the bath from the surround- 
ings. 

Numerical result with inertial term added 

In the following we are going to present some results ob- 
tained from the direct numerical simulation of the full 
model containing inertial terms and we will see that, the 
basic observation of the over-damped (simple) model is 
also present there. We will see that, the velocity of the 
CM scales as squire root of the force constant a and ksT 
and picks up the write direction as is evident from the 
discussion of the symmetry breaking. The absolute value 
of the velocity is somewhat smaller than that obtained 
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FIG. 6. Plot of < V > vs a. 



FIG. 7. Plot of < V > against a. 



from the over-damped limit for obvious reasons of having 
velocity fluctuations considered in the full inertial model 
in the presence of the masses. However, an analytical 
solution of the full model with inertial terms, till now, 
has not been possible to do because of the strong nonlin- 
earity introduced by the collision potentials. The model 
is the following 
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While simulating the above model, we will exchange the 
velocity and position of particles at each overshooting 
considering them to be of the same mass. This is the 
implementation of the hard collision happening in the 
absence of the damping as has been explained in details 
in the context of the over-damped model. For the full 
model, we are also exchanging velocity on top of doing 
the same for the position. This is a very simple way of 
implementing the basic requirement of our prototypical 
model of having non-simultaneous working of the damp- 
ing (slow) and the collision (hard) interactions in the nu- 
merical simulation which actually proves technically dif- 
ficult in the analytical handling of the full inertial model. 

In Fig. 6a and 6b we have plotted the average veloc- 
ity of the CM against a. Fig. 6a shows plots for the 
value of k B T = 0.5 (lower one) and k B T = 0.005 (up- 
per one), whereas, Fig. 6b shows plots of the same for 
k B T = 0.00005 (lower) and k B T = 0.0000005 (upper). 
Its clear from these figures that the system is showing 
similar results as the over-damped system but the ab- 
solute value of the CM velocity is generally smaller by 
2 orders of magnitude than the over-damped one. The 
masses taken here are unity as is obvious from the model 
mentioned above. The simulations have been done with a 
time step At = 0.01 for 10 8 steps for every set of param- 




FIG. 8. Scaling of | < V > | with a. 



eters and then averaging over 50 such runs each starting 
from the same initial conditions. If we take the mag- 
nitude of the velocity and plot the graphs on log scales 
against a (Fig. 7) and k B T (Fig. 8) we immediately sec 
that the velocity scales as \fa and \/k B T as is predicted 
by the over-damped model which confirms the validity 
of our exact over-damped analysis handling two different 
time scales. Actually, the exponent of a varies between 
0.5 and 0.6 but we think a better averaging would fix it 
close to 0.5. The scaling with a remaining unchanged is 
important in view of the fact that its the harmonic in- 
teraction which basically drives the system. The scaling 
with y/k B T is important in view of the fact that the sym- 
metry broken ID system being bound to move in a single 
direction should also respect equipartition in equilibrium 
which is a symmetry of the equilibrium itself. The de- 
viation in the amplitude of the velocity from the over- 
damped expression cannot be fixed without being able 
to do the analytic work on this model. Anyway, having 
a basic confirmation of important scalings from the nu- 
merics of the full model lends support to our simplified 
analysis with the over-damed model. 

Its important to note that, standard equilibrium 
statistical physics is actually in contradiction with 
Newton's first law if it totally denies the existence of 
any uniform motion under the force balanced situation. 
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FIG. 9. Scaling of | < V > | with k B T. 



A mechanical equilibrium in general is a force balanced 
situation and some velocity gained during the initial 
transient to equilibrium should be maintained if under 
certain conditions equilibrium renders the force balance. 
Our Eq.6 is the statement of such a force balance where 
the damping of the velocity of the CM is equal to 
the average force acting on the system due to broken 
symmetry. Its kind of a self-consistency condition that 
selects the average velocity of the CM. Just because 
such a situation had previously not been conceived of, it 
cannot be enough reason to think that there cannot exist 
a situation under equilibrium which would be consistent 
with Newton's first law specially when the Hamiltonian 
dynamics is the basis of equilibrium statistics. The 
problem so far has been that the results of particle 
dynamics imposed to everything. We see that, the 
moment we start considering the structure within the 
particle, there are more scopes of having new results. 
Moreover, equilibrium being defined by a whole lot of 
stationary quantities, one should keep in mind that 
those stationarity are average stationarity. Thus, there 
is scope of having other structures withing the average 
time period over which this stationarity is defined. One 
such possibility has been discussed in the relation with 
inelastic collisions and possible relatively long lived 
fluctuations bringing in so far unknown results. In 
general, it appears interesting to look at the effect of 
these internal geometric and functional structures in 
the fluctuations within equilibrium, specially when, we 
are able to do experiments at this scale these days. 
It probably also would provide means of studying 
non-equilibrium and equilibrium on the same footing 
because when fluctuations are long-lived they actually 
enable us seeing non-equilibriums within equilibrium. 
Discussion 

Having understood the physics at the origin of sus- 
tained uniform motion in equilibrium let us emphasize on 
the fact that such a directed motion is perfectly consis- 
tent with the Kelvin statement of the second law of ther- 



modynamics. Since, we are considering a closed system 
of the bath and the dimer, no energy is being extracted 
from the system and the situation in no way violates en- 
ergy conservation. Any extraction of energy from such 
a system would start a transient where the bath would 
start losing temperature because no bath has infinite heat 
capacity and the system is no longer in equilibrium. Such 
a transient would persist up to the point where the ve- 
locity of the system (being proportional to the square 
root of the temperature) is not enough to drive it against 
a barrier and in the resulting equilibrium situation the 
dimer would come to rest. So, in no way there could be 
perpetual extraction of energy making the dimer move 
against a global barrier. Our model does not include 
such a situation, its an isolated bath dimer system not 
having any other external drives, barriers and dissipative 
interactions. 

Let us also note that for a uniformly moving system 
of particles in equilibrium the Detailed Balance would 
also hold good due to the fact that homogeneity of space 
makes interactions between particles function of the sepa- 
rations between them. So, no contribution to the internal 
energy of the system is there due to CM motion. Equilib- 
rium probability distribution, that ensures detailed bal- 
ance, being a function of the interactions between parti- 
cles would remain unaltered. This fact is explicit in the 
form of P(Z) we have got. Moreover, the heat bath does 
not set any reference frame to the dimer means the prob- 
ability of the CM to be found anywhere in the bath is 
a constant. This, in other words, is the homogeneity of 
space for the dimer. The homogeneity of space is also im- 
portant to understand that no entropy production would 
be there for such a system under uniform motion in equi- 
librium. The rate of entropy (S(t))) production is given 
by the expression 

Where P(x,t) is the probability to have the configura- 
tion x at time t. Now, the configuration being defined 
solely on the basis of the separation between particles and 
keeping homogeneity of space in mind (as our P(Z)), the 
P{x, t) would never be an explicit function of time in 
equilibrium i.e. P(x,t) = P(x). Thus, the partial time 
derivative of the configuration probability being zero the 
entropy is a constant which is the most important sig- 
nature of equilibrium from which basically follows the 
detailed balance etc. Physically one can understand this 
constancy of the entropy from the constancy of the avail- 
able phase space to the system. The conservation of the 
allowed momentum space to the system remains there 
in the presence of a constant drift of the CM. Actually 
the Maxwellian velocity distribution is centered around 
the CM velocity in such a situation and not around zero. 
The available volume to the system also remains con- 
stant since, due to uniform motion of the CM, as much 
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new volume is getting available to the system as is lost 
at every infinitesimal time interval. This also keeps the 
available volume to the bath degrees of freedom constant. 

We first understand why the Smoluchowski-Feynman 
ratchet does not show directed transport in equilibrium 
whereas our system does. The reason of Smoluchowski- 
Feynman ratchet not showing directed transport lies in 
the fact that it is not at all a symmetry broken system 
when the temperatures of the two heat baths are the 
same i.e. in equilibrium. This fact has beautifully been 
demonstrated by Feynman showing that the probability 
of the pawl going up is the same irrespective of which 
bath it is driven by when the temperature of the baths 
are the same and that removes the broken symmetry [j| . 
But, our system is a broken symmetric system perfectly 
consistent with the conditions of equilibrium and that 
is why it shows directed transport at constant tempera- 
ture. So, directed transport in the absence of a driving, 
being in contact with a heat-bath, is possible. It does not 
automatically mean available energy and one should not 
confuse our present results with perpetual motion ma- 
chines. Note that there exists example of uniform motion 
in equilibrium namely in superfluid helium and supercon- 
ductors where a gap in the spectrum keeps the system 
away from excitations. So, one can reasonably under- 
stand that had the uniform motion and thermodynamic 
equilibrium been generically inconsistent because of de- 
tailed balance, entropy production etc., we would have 
not got such states. So, there exists no generic reason as 
to why one cannot expect such a situation in mechanical 
systems as well. Rather, our model shows how the de- 
tailed balance ensures broken symmetry, hence, directed 
motion keeping the energy balance intact. 

Le us have a look at very general dynamical system and 
its corresponding Fokker-Planck equation to have some 
clue as to how a directed motion is not at odds with the 
the demand of equilibrium. This part of the discussion 
would also help us understand what it actually means 
that the bath does not set any reference over space to 
the system, consider the most general form of dynamics 
of a particle as 



dx{t) 

dt 
dv(t) 

dt 



v(t) 

Tv(t) - 



dV{x) 
dx 



•/(*), 



(10) 



where in the above model, V(x) is the potential account- 
ing for all the forces other than the stochastic force f{t) 
and the damping where T is the damping constant. The 
Fokker-Planck equation for such a system can be imme- 
diately written (in ID) as 

dP{x,v,t) dP(x,v,t) 1 8V(x) 8P(x,v,t) 



dt 
m dv 



vP(x, v, t) 



dx to dx 

k B T dP(x,v,t) 

777 dv 



dv 



(11) 



This Fokker-Planck equation has a solution of the form 
P(x,v,t) = Nexp(-V(x)/k B T - mv 2 /2k B T) where N 
is the normalization constant. Had we done a Galilean 
transformation to our system Eq.10, its easy to see that 
we would have got a similar Fokker-Planck (due to v 
appearing linearly in Eq.10) with v replaced by v + C 
(where C is the constant velocity) and the solution of the 
Fokker-Planck in that case would have been P(x,v,t) = 
Nexp(~V(x)/k B T ~m(v + C) 2 /2k B T). So, the prob- 
ability distribution has a uniform drift. Important to 
note that this uniform drift is actually through the bath 
although, of course, one would see the minimum of the 
global field V(x) to move with the same velocity and the 
particle equilibrating at this minimum. But, this motion 
through the bath is something important to note which 
appears because of its homogeneity and of course not 
noticeable because of the same homogeneity. Neverthe- 
less, this proves a point that under the above mentioned 
consideration of the dynamical presence of the heat bath 
with a damping and fluctuating force the rest of the sys- 
tem can actually move through it without experiencing 
any additional damping because of the baths homogene- 
ity. This also indicates to the fact that such uniform 
motions can also be possibly present when the field is 
not fixed in space of the observer, but, can actually move 
with the system i.e. when the field V(x) is due to in- 
ternal interactions between a system of particles and has 
a minimum for the system to equilibrate in. Note that, 
the structure of the Fokker-Planck in such a case should 
remain the same with some additional terms on the right 
hand side due to the fact that internal field can now 
introduce some correlations between the velocities and 
a selection of the average velocity of the system should 
come from these additional terms. For our system with 
a hardcore collision i.e. not a smoothly varying confining 
field V(x) it appears difficult to derive the related Fokker- 
Planck incorporating the velocity exchange at collision. 
Also note that, with inelastic collision and smoother col- 
lisions also its difficult to get the Fokker-Planck due to 
the extra dissipation present. However, We at present are 
doing work on it hopefully to get useful results which we 
can publish later. But, our present analysis being simple 
actually captures such a possible situation which is very 
much plausible in quite an easy way. 

The most important part in the construction of our 
model is the consideration of two time scales where at 
the slow scale the symmetry is broken and the fast scale 
is symmetric. This two time scale scenario can be con- 
sidered as a generic mechanism for construction of such 
models. Important to note that the Langevin picture de- 
veloped for a Brownian particle is based on integrating 
out the faster degrees of freedoms of the bath. But, in 
the case of two Brownian particles bound to each other 
where they frequently collide with each other and the 
collision time is much small compared to the Langevin 
time scale - that part of the fast dynamics has not been 
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integrated out. It seems a simple but extremely effec- 
tive technique to take into account the faster dynam- 
ics through the boundary conditions. The correctness of 
this technique ultimately has to be verified by experi- 
mental results. Note that, treating simultaneously two 
such widely different time scales and forces where the 
dynamics essentially becomes nonlinear is extremely dif- 
ficult if not impossible. Appropriate boundary conditions 
based on experimental observations is one possible way 
out and here we are proposing one which definitely has to 
be checked by experiments. The source of our conviction 
at the present case is (a) the excellent agreements with 
the numerical results (b) collision does not feed energy to 
the system so that < Z > can go on increasing for ever. 

In the present work we have argued that second law 
and detailed balance hold good. We have shown that the 
entropy remains constant as it should be in equilibrium. 
The only difference with the so far existing knowledge of 
equilibrium is that the time reversal symmetry present in 
stationary equilibrium is not there for uniformly moving 
systems. We know consequences of second law that fixes 
an arrow of time goes against time reversal symmetry. 
Its also important to note that the master equation for 
the probability distribution dynamics does not have time 
reversal symmetry. Then, why only a restricted class of 
solutions of the master equation with time reversal sym- 
metry has to be considered as equilibrium solution par- 
ticularly when uniformly moving ones are also entropy 
conserving? Moreover, a uniform motion is a mechanical 
equilibrium situation. It does not contribute to the tem- 
perature of the system because temperature is a measure 
of the width of velocity distribution. So, thermal equi- 
librium does not get perturbed by uniform motion of the 
whole system. There can indeed be stationary solutions 
of the master equation on uniformly moving frames which 
can always equally qualify as equilibrium solutions of the 
system having the time reversal asymmetry of the basic 
dynamics (master equation). An important point to be 
noted is that, within the scope of our present treatment, 
we never know what happens when the hydrodynamics 
of the heat bath becomes important because Langevin 
dynamics does not take that into account. Our analy- 
sis, as many others, is based on the assumption that in 
over-damped regime the effect of hydrodynamics is neg- 
ligible. Also note that, under non-equilibrium conditions 
when we get directed transport (there are numerous pa- 
pers on such driven diffusive systems) we take the driving 
to have no directional bias (that's what makes the result 
nontrivial) and still we get directed motion of the object 
under consideration in contact with the heat bath. So, a 
conservation of linear momentum would require an over- 
all drift of the bath degrees of freedom in the direction 
opposite to that of the transporting system. If this effect 
is to be important to the extent it prevents the directed 
transport we have seen under equilibrium conditions, so 
it would do in the non-equilibrium conditions as well. 



To conclude, we would like to say that, with the help 
of a somewhat idealized model we probed the problem 
of having symmetry breaking directed motion in equilib- 
rium because the model is exactly solvable. The phys- 
ical understanding we gain from this analysis leads us 
to realize possible conditions under which such a motion 
can be physically realizable. We understand that dissi- 
pation plays a crucial role in such processes and based on 
which we can actually relax the stringent conditions of 
our idealized model. A possible schematic description of 
an experiment would be the following. If it is possible to 
implant a charge on one of the constituent particles of the 
dimer having the effective attractive interaction between 
the particles of the dimer by other means (e.g. effective 
hydrophobic interaction) then one can use a magnetic 
field to have such systems rotate on a circle while hav- 
ing some average velocity of the CM. Remember that 
the Lorenz force would not do any work on the system. 
Most probably one would have to give the system some 
initial velocity just to break the isotropy in the begin- 
ning and then the system should eventually settle to its 
average speed under the given conditions of the bath. 
We strongly feel that this is an area worth investigating 
because it can change the way we understand statistical 
physics at the present times. 
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